Kondo effect and singlet-triplet splitting in coupled quantum dots in a magnetic field 
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We study two tunnel-coupled quantum dots each with a spin 1/2 and attached to leads in the 
Coulomb blockade regime. We study the interplay between Kondo correlations and the singlet- 
triplet exchange splitting K between the two spins. We calculate the cotunneling current with 
elastic and inelastic contributions and its renormalization due to Kondo correlations, away and at 
the degeneracy point K = 0. We show that these Kondo correlations induce pronounced peaks in the 
conductance as function of magnetic field B, inter-dot coupling to, and temperature. Moreover, the 
long-range part of the Coulomb interaction becomes visibile due to Kondo correlations resulting in 
an additional peak in the conductance vs to with a strong B-field dependence. These conductance 
peaks thus provide direct experimental access to K, and thus to a crucial control parameter for 
spin-based qubits and entanglement. 
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Semiconductor quantum dots have attracted much 
interest over the years due to their rich and repro- 
ducible transport properties in the Coulomb blockade 
(CB) regime, where the number of electrons on the dot 
is quantized due to charging effects |Q. More recently, 
attention has focused on the spin of the electron in such 
nanostructures, with the spin introducing new correla- 
tion effects such as Kondo behavior ^~^], and its in- 
terplay with spin exchange interaction in single fi] [n]] 
and double dots Jll],[l^]. On the other hand, it has been 
pointed out that the spin, confined to a quantum dot 
or atom, is a suitable candidate for quantum informa- 
tion processi ng j]l3| |, due to unusually long decoherence 
time of spin jl4| . A crucial element in such spin-based 
quantum computing schemes is the Heisenberg exchange 
interaction K (singlet-triplet splitting) between spins of 
adjacent dots, being controlled via the interdot tunnel- 
ing [Q. Thus, the primary goal is to achieve control 
over K which then allows one to generate deterministic 
entanglement of spins. Using a Hund-Mulliken (HM) ap- 
proach to describe a realistic double dot system (DD) 
it has been shown that K is very sensitive to long 
range Coulomb interaction as well as to magnetic fields 
by which a singlet-triplet crossing can be tuned. Moti- 
vated by this we study here transport and Kondo effects 
in such a realistic DD system within the HM approach, 
thereby going beyond short-range on-site models used so 
far to describe Kondo effects. In particular, we calcu- 
late the current through the DD via a Schrieffer- Wolff 
transformation and via a systematic cotunneling calcu- 
lation including elastic and inelastic contributions. Us- 
ing a perturbative renormalization group (RG) approach 
we show that the conductance in the cotunneling regime 
shows pronounced peaks induced by Kondo correlations 
and long range Coulomb interactions as function of tem- 
perature, inter-dot coupling, magnetic fields, and bias. 
Such Kondo enhanced peaks in the conductance thus pro- 
vide direct experimental access to singlet /triplet states 
and their exchange splitting X-the quantities of crucial 
importance for spin-based qubit schemes. 



We consider a DD system consisting of two lateral 
quantum dots tunnel-coupled to metallic leads, in the 
presence of a perpendicular magnetic field B, see Fig. 1. 
At low temperatures T, the conductance G of the DD 
as a function of the gate voltage V g shows sharp dou- 
blets of sequential tunneling peaks separated by CB val- 
leys (cotunneling regime) p™. In the middle between 
such peaks, the number of electrons in the DD is even, 
M = 2N (assuming similar dots). We assume that N — 1 
electrons of each dot form a closed shell [with N—l even] 
and thus ignore them. The remaining two electrons in the 
DD can be described by the Hamiltonian 



i=l,2 



A(r,) +W{r t 



+ C, (1) 



where C = e 2 /«|ri — ra| is the Coulomb interaction, with 
charge e and dielectric constant n (= 13.1 for GaAs). As 
usual in GaAs, the Zccman interaction will be neglected. 
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FIG. 1. Double-dot system containing two electrons and 
being coupled in series to two metallic leads at chemical 
potentials firt and fj,L with bias Afi = — fiR. The 
electron spins Sl, Sr interact via the exchange interaction 
K = E t — E$, where Et,s is the triplet /singlet energy. 

However, the orbital effect of the magnetic field B = 
V x A is important, it allows us to tune the singlet-triplet 
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transition in the isolated DD, provided, however, there is 
long-range Coulomb interaction . Within the low en- 
ergy sector, a realistic confinement potential for the DD 
is given by g| W(r) = {mm 2 /2 )[{x 2 - a 2 ) 2 /4a 2 + y 2 ], 
which separates at a > yjTi/mu>o into two harmonic wells 
of frequency u>q, located at x = ±a, with, typically in 
GaAs, hujQ ~ 3 meV JjJ. Within the HM approach we 
express the two-particle states of the DD jL5| in terms of 
symmetric (for spin singlet) and antisymmetric (for spin 
triplet) combinations of the Fock-Darwin states, given by 
tp±a(x,y) = c xp[-((s T a) 2 + j/ 2 )/2 A 2 T iya /2l 2 ]/Xy^, 
where I = y/hc/\e\B and A = Wh/mu}, and ui = 
+ lo 2 l , where cj^ = |e|£?/2mc. We introduce the 
operators d± a {d±^ a ), which create (annihilate) the DD 
states V±,<j = Xa(<^-a ± V+ a )/v / 2(l ± S), where x<r is 
the spinor and 5 = (<y2±a|</?=Fo) the overlap integral. The 
singlet and triplet states then become, 



|00> = 7/r+F (< ^ +; ~ Hl ^ m ' 
|ii) = dL T 4 T |o) , |i-i) = dLXil ) » 



(2) 



|io) = -^=(dL T 4 i + dL i < T )|o), 



t jt 



where |0) is the DD ground state with M = 2N — 2, and 



1 + 



(3) 



where t R = to + tc is the extended inter-dot tunnel- 
ing amplitude with to its bare value (C = 0) and ic(to) 
the part that gets renormalized by the Coulomb interac- 
tion C. Similarly, Uh is the extended on-site Coulomb 
repulsion which decreases with decreasing a Note 
that for to — we have 0=1, while 4> < 1 occurs 
due to double occupancies in the dots, and 4> ~~ * for 
C — > 0. The DD is connected to two Fermi liquid leads 
a = L, R (Fig. 1), described by H = J2 a ka £ kC f aka c aka , 
where c^ aka creates an electron-state with momentum k 
and spin a in lead a. The tunnel coupling from lead a 
to the nearest dot is parametrized by the amplitude t a , 
and thus the amplitude to tunnel from lead a into the 
DD states ip na is then given by t^,± = tt,/y2(l ± 5), 
t R< ± = ±t R /y / 2(l ± S). We use the tunneling Hamil- 

tonian H T = I], iQfco .(ia n c^ fe(T d, i(T + h.c.) to perform a 
Schrieffer- Wolff transformation [jl7| , and arrive at the co- 
tunneling part H co i = ifr[(*0 + — Lq)^ 1 Ht], with Liou- 
villean LqA= [Ho, A], where Hq = Hd + Hi. 

Next we concentrate on the states (^) and project out 
the higher lying energy levels by means of the mapping @] 
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(4) 



where S± = S L ± S R , T = S L x S R , and <j>± = y/2(l ± 
4>)/(l + <fi 2 )- The spin 1/2 operators Sl.r represent the 
two electron spins on the DD We arrive at the 

effective Hamiltonian 

H = H+KS L -S R + AH, (5) 

+/ Q + Q -s Qa - S_ + 2iI- a ,s aa ,-T) , (6) 

where s QCe / = £ fcfcW c' aka (a aa > /2)c a i k ' a ', and p aa / = 
Efcfc'tr c Lfc CT c oi'fe'CT • In @ the bare constants J, V, I ± are 
defined at the energy cutoff Dq ~ Tiwo, 



2E r 



-7" 



T, 



(7) 



with matrix elements , = f * 

W^a' = £',-*«■+ ± &,+*«,-, and £ c = 2E+E_/(E+ + 
E-). Here, £± = £(M ± 1) - E[M) T- M is the CB ad- 
dition/extraction energy, with E(M) being the energy of 
the DD with M electrons, and fi = {fiL + We 
calculate E(M) for M =1,2,3 within the HM method. 

Attaching leads to the DD results in a shift of the DD 
spectrum such that K = Kq + 5Ksw > where 



Ko = v - 



Uh 
2 



yu H +m H 



(8) 



is the exchange interaction of the isolated DD, with 
v accounting for long range Coulomb effects and 
SK SW = 2vD (E c /E^)TtV < stems from the 
Schrieffer- Wolff transformation. We next derive the RG 
equations for the coupling constants in ^) by applying 
the "poor man's" scaling approach p8[ . In matrix form 
the RG equations read, 



j=J 2 + {I+)' 



- (I-f 



V = 2[I-,I^ 



(9) 



where the dot denotes d/d(v£). Here v is the density of 
states per spin in the leads, and £ = \n(Do/D), with 
D being the scaled cutoff. Eqs. ^ are valid for D S> 
K, A/x, T (for T see below). The set of 16 equations (||) 
can be reduced to 9 equations for our case of t an being 
real jl9). Solving ^ numerically we find the character- 
istic energy scale of the problem, To = Dq exp(— 7/1/ Jo), 
where Jo = (t 2 + t 2 R )/Ec, and 7 < 0.5 is a non-universal 
number, determined by the ratios of the coupling con- 
stants (^) to Jo, and thus depends on the internal features 
of the DD. We note that at t R = we find 7 = (l-S" 2 )/4, 
i.e. independent of the Coulomb interaction. Also, in the 
limit of separated dots, to — 0, Eqs. (||) decouple into two 
RG equations, determining the Kondo temperatures for 
two 1/2 spins @, T« K{L R) = D exp(~E c /U 2 L R ). 

We derived Eqs. @ assuming a singlet-triplet degen- 
eracy (STD) of the DD, i.e. K = 0. However, the RG 
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procedure generates terms that also renormalize the ex- 
change K between the DD spins. Indeed, we find 

K = 21n(2)L>t/ 2 Tr{V 2 + (I + ) 2 - (rf - J 2 } . (10) 

The RG flow (|Io|) resembles the generation of the RKKY 
interaction in the two-impurity Kondo model p0| . It 
follows from ([l0]), (||) that, at low T ~ Tq, the renormal- 
ized K = K(B, D) strongly depends on both D and B 
at the STD point with K(B^ 1 D) = 0, see right inset 
of Fig. 3. This suggests that, (i) B* shifts towards the 
triplet side as T is lowered down to To, (ii) the energy 
scale of the problem strongly depends on Kq , which pre- 
sumably implies a rather sharp crossover between a spin 
1 Kondo regime on the triplet side and a locked singlet 
of Si, Sr (i.e. no Kondo effect) on the singlet side. 

We perform a non-equilibrium calculation of the cur- 
rent through the DD up to the second order in the pertur- 
bation (||) . The current consists of an elastic and inelastic 
component, T K = l el + l lne i, 



lel = 

n 



J 2 lr(S 



V, 



LR 



Am, (ii) 



l m ei = —™ 2 (lZ R ) 2 g(PK)f{pK, /3AM) , (12) 

where (3 = 1/T (with Boltzman constant set to one), 
g(u) = (e u + l)/(e"/3 + 1) accounts for the degeneracy 
of the excited level, (Si) = 6/(e /3K + 3), and 



ii tanh(it/2) sinh(w) + v[l — cosh(w)] 
cosh(w) — cosh(v) 



(13) 



Eqs. (|i~l"|), (|l2|), together with the RG equations (||), de- 
scribe the renormalization of the current 1, i. e. I oc 
1k(^(Do/1))- However, not all terms were retained in 
(|J), but only those which become renormalized. To give 
a complete description of 1 at temperatures Tr- <C T < 
Dq, we perform a systematic calculation of the cotun- 
neling current using the technique developed in pjj . In 
the cotunneling regime, including level shifts but neglect- 
ing heating effects pjj, we match the two currents and 
obtain for the total current (after lengthy calculation) 
1(C) = 1 K (£) + 61(0) , with 



61 



V? 



2VlrJlb, 



LR 



<s 2 + ) 



(2-(S 2 + )) 



J 2 



LR 



E_ 



E- + Ej. 



Ali 



(14) 



In the cotunneling regime, the differential conductance 
G = edl/dAfj, shows a step at A/i = ±K. for T < \K\. 
This step is due to the inelastic current ( |12| ) contributing 
for | A/i | > \K\. The step height of G is different for 
the ground state being a singlet or triplet; the ratio of 
heights is given by g(K/T). As T is lowered, Kondo 



correlations are expected to develop at A/i = ±K, with a 
Kondo peak arising at the step of G. Kondo correlations 
at finite bias are not the subject of the present paper, 
but we would like to point out that the step in G can 
be used for determining the ground state of the DD and 
measuring K. 




FIG. 2. Linear conductance G(to) at different B's. Dotted 
lines: cotunneling contributions. Dot-dashed line: G vs to 
at B = OAT without the long range part of the Coulomb 
interaction, i.e. for v = in (H). For definiteness, we keep 
E+ = E- (middle of CB valley) by adjusting V g when varying 
to, and set if, — tn in Figs. 2-4. Inset: Comparison of full 
(solid line) with chain (dashed line) model Hamiltonian. 

We turn now to a discussion of the linear conductance, 
G = G(Afi = 0). First, we consider the case B — 0. At 
small to the RG growth of K is weak (due to near cancel- 
lation of the trace terms in ( |i"o| ) ) , and a limit is reached 
where each spin is strongly coupled with a lead electron. 
The Kondo effect of the DD, in this case, consists of two 
independent spin 1/2 Kondo effects for each of the dots 
separately. At large to, the RG correction to K can be 
neglected because of a large value of K . The two dot 
spins are locked into a singlet state in this case, and the 
lead electrons feel only the potential scattering, slightly 
enhanced by the RG flow (^) which terminates at D = K. 
At intermediate to, such that K ~ Tk, the exchange in- 
teraction (between the dot spins) and the spin 1 /2 Kondo 
effect compete, and a crossover between the two regimes 
occurs. At this crossover, the renormalization of K is 
comparable with the Kondo temperature Tk\ each of the 
dot spins couples to both leads. The conductance versus 
to shows a peak at to ~ ^/UhTk (see Fig. 2). A second 
peak at larger to emerges due to the interplay between 
Kondo correlations and long range Coulomb interaction, 
see dot-dashed line in Fig. 2. This striking sensitivity 
of the second peak on the long range Coulomb interac- 
tion provides a way of studying the screening effects in 
coupled quantum dots. As a check we have taken the 
Hubbard limit (chain Hamiltonian) of our model, and 
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FIG. 3. Linear GvsB around K = 0. The kinks in the dot- 
ted-line regions are an artifact of our two stage RG procedure 
in the crossover regime and will be smoothened in an exact 
treatment. The stars denote -B* with K(B,,D — Do) = 0. 
Left inset: Cotunneling conductance vs B at small coupling to 
the leads where Kondo correlations become negligible. Right 
inset: Renormalized K vs B, showing a sharp singlet-triplet 
transition at low T(~ To). 

found good agreement with the exact NRG calculations 
of Ref . jl2| . In inset of Fig. 2 we compare the two models 
to illustrate the importance of the long range Coulomb 
interaction. For small B (with K > 0), we find that the 
second peak in G(t ) is very sensitive to B, see Fig. 2. 
0.6 




T/D 

FIG. 4. Linear G vs T for different to. Insets: Compari- 
son of cotunneling (dashed line) and renormalized (solid line) 
calculations for different couplings to the leads. 

At larger B, a singlet-triplet transition occurs jul . 
Around this point K = 0, the RG flow ([|) terminates at 
D = T > | if |, which leads to a striking RG enhancement 
of the peak in the cotunneling conductance, see Fig. 3 
and left inset. In the singlet regime, for K > T, the RG 
flow terminates at D = K, and there is no Kondo limit. 
Physically, this means that the DD spins are locked into 



a singlet with exchange energy K > 0. In the triplet 
regime, for —K > T, we apply a two stage RG procedure 
p2"t and find a i\-dependent energy scale for the spin 1 
Kondo effect, very much like as in single dots [|]||. At 
low T, the renormalized K shows a striking dependence 
on B around K = 0, see right inset of Fig. 3. Our per- 
turbative RG thus suggests a sharp transition between 
the two regimes, with a rapidly vanishing spin 1 Kondo 
temperature when approaching the singlet side. Finally, 
in Fig. 4 we plot the conductance vs T. We see that again 
the Kondo correlations lead to a pronounced peak in G 
which occurs at T = K , thus providing a further way to 
access the exchange interaction K . 
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